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ABSTRACT 

The Nambu - Jona-Lasinio model in its SU(2) and SU(3) versions with scalar 
and pseudoscalar coupling are applied to baryons. The parameters of the model 
are fixed in the meson sector. The baryons arise as a soliton of three valence quarks 
coupled to the Dirac sea (quark-antiquarks pairs). Within the SU(2) version the 
nucleon static properties as well as some form factors, namely, the electric and 
axial ones are described quite successfully. The nucleon-delta splitting comes out 
reasonable. In medium the nucleonic soliton gets less stable - the mass is reduced 
whereas the radius increases. At some critical medium density there is a clear 
derealization of the soliton: The nucleon does not exsist anymore as a soliton. In 
SU(3) version the strangness carrying baryons are described as SU(3)-rotational 
excitations of the SU(2)-soliton embedded in the SU(3)-sector. The mass splittings 
between the octet and decuplet as well as within the multiplets are peproduced 
not only in a correct oder but also in a good agreement with the experimantal 
values. 
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1. Introduction 



In the last years because of the tremenduos difficulties in the nonperturbative 
(strong coupling) regime of QCD there is an increasing interest in the effective 
chiral models respecting the chiral symmetry breaking. Among these models the 
Nambu - Jona-Lasinio model applied to quarks is the only one which allows to 
combine in a natural way the two extreme opposite pictures - the pure valence one 
advocated by the constituent quark model and that used in the Skyrme model. 
The model makes use of the experimentally supported "standard" picture of the 
nucleon as a bound state of three valence quarks coupled to the mesons arising as 
quark-antiquarks pairs from the Dirac sea. Similar to the other effective models the 
NJL model posseses two shortcomings. On the one hand in the model the quark 
confinement is not implemented and one can only hope that the confinement is not 
of great importanse for the very low-energy baryon properties. On the other hand 
the model is not renormalizable. It means that a finite cut-off is needed to make 
the theory finite. The cut-off introduces into the theory a physical scale which 
should be fixed by reproducing some physical quantity. 

2. NJL model with SU(2)-fiavour 

We start with the simplest lagrangean of the NJL model with scalar and pseu- 
doscalar quark-quark couplings 1 : 



Here * describes a quark field with SU (2) -flavour (u and d) quarks, N c = 3 colours 
and the average current quark mass mo = (m u + m^)/2. Introducing auxiliary 
sigma and pion fields by a = —g^^/X 2, and n = — g^i'j^r^ / X 2, and assuming 
them to be classical (zero boson and one fermi loop approximation), one can use 
path integral techniques in euclidean metric to evaluate the corresponding effective 
action as a sum of a quark part 



S q eff = Sp\og — + h-fji) with h = — + /3g(a + i^n.r) (2.2) 




(2.1) 
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with n being the chemical potential, and a meson one 




(2.3) 



2 



In the meson part using the PCAC the current mass tuq is eliminated in favor of 
the pion mass and the pion decay constant f n . The coupling constant G is 
related the new one A by G = # 2 /A 2 where the additional coupling constant g is 
introduced for convinience. The latter will be fixed later to the physical coupling 
constant in order to identify the auxiliary pion and sigma fields with the physical 
ones in the meson sector. Only one part of the effective action, S^^{ji = 0), 
coming from the Dirac sea is divergent. Using the proper time scheme 2 we make 
it finite. Subtracting the vacuum contribution one can express the energy (see for 
more details refs. 3 ' 4 ) as a sum of a quark part 



and the meson part which up to a trivial factor is given by eq.(2.3). The energies 
e a are the eigenvalues of the hamiltonian (2.2) and correspond to the vacuum 
solution. 

We fix the parmeters of the model, namely the coupling constant A 2 and the 
cut-off A, in the meson sector of the vacuum by reproducing the experimental 
values of the pion decay constant f n and the pion mass m n . This procedure 
originates from Eguchi 5 and is described in detail in ref. 3 . We use the stationary 
meson field configuration < tx > v = tv v = and < a > v ^ where the the latter is 
given by the non-trivial solution of the well-known gap equation. It generates the 
constituent quark mass M = g < a > v = gf n . Actually we use the gap equation 
to express the coupling constant A 2 as a function of the constituent mass M. For 
the meson dynamics one can use the meson propagators^ and to fix the g to the 
physical pion coupling constant in order to identify the auxiliary a and n with the 
physical ones. The next step is to reproduce the physical value of the pion decay 
constant. The latter leads to a relation between the M and A. However, in the 
present case one can also use the gradient expansion which means that the coupling 
constant is fixed to yield the correct form of the meson kinetic energy term and 



E«= £ e a + 



0<e Q </x 




(2.4) 



with the proper-time regularization function 




(2.5) 
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the pion mass is related to the second derivative of the energy at the stationary 
point. It should be also noticed that in this limit (q 2 = 0) the pion decay constant 
is reproduced trivially for the physical pion fields. All those conditions together 
leave the vacuum value of the constituent mass M as the only free parameter. In 
principle, the empirical values 8 of the quark condensate and the quark bare mass 
can be used to fix M but actually they still leave a broad range for M. 

For the solitonic sector we solve the Dirac equation (eq.(2.3)) with the equa- 
tions of motion of the meson fields: 

^ = >« E^1/2(^ A )- E Ma\-^. (2.6) 
A I a 0<e a <^ J 

71 = X2 Nc { E<^75(-r-f)<M#l/2( e a> A ) - E 0a*75(™)0a > . (2.7) 
A I a 0<e a <v J 

The meson fields are assumed to be in a hedgehog form and are restricted on the 
chiral circle a 2 + n 2 = / 2 . We use a numerical self-consistent iterative procedure 3 
based on a method proposed by Ripka and Kahana 9 . 

As a next step the classical soliton should be quantized. To that end one 
can make use of the rotational symmetry of the solution introducing isorotating 
hedgehog meson fields 

U(x) -> R{t)U{x)R\t), (2.8) 

where R(t) is a unitary SU(2) matrix. Following the usual canonical quantization 
procedure 10 one can relate the angular velocities to the spin (isospin) operators 
by 

en l (n a ) -> -u i {T a ). (2.9) 

Thus one can assign proper spin and isospin quantum numbers to the soliton and 
evaluate the nucleon properties like energies, radii, form factors, spin and isospin 
contents, etc. 

For the energy an expansion up to the first non-vanishing order with respect 
to angular velocity gives the well-known rotational energy correction to the clas- 
sical energy from which one should subtract the corresponding translational and 
rotational spurious zero-point energies 11 

f'-f 1 J ( / + 1 ) <t2> <p2> r2 im 

E - E «* + —20 29 29— (2 - 10) 
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Here is the moment of inertia given as a sum of valence and sea part . The 
latter is reguralized in the proper time scheme. The results of the energy of the 
nucleon and the delta are shown in Fig. 1 and some values of the moment of 
inertia for different constituent masses are presented in Tab. 1. As can be seen the 
experimental nucleon-delta mass splitting is reproduced well for M around 420 - 
450 MeV. The sea contribution is less than 20 %. 

The nucleon electric form factors are given by 



where \N(p)) is the quantized nucleon state and jg 171 is zero-component of the 
electromagnetic current operator defined as 



To that end an external electromagnetic field is coupled. Similar to the moment 
of inertia the form factors incude also a valence and sea contribution. It should 
be noted that the isoscalar form factor does not need any regularization which is 
not the case of the isovector one. The proton and neutron electric form factors 

i o 

calculated in the non-relativistic approximation q <C Mjy- are shown in Fig. 2 
and Fig. 3, respectively. The corresponding neutron charge distribution presented 
in Fig. 5 has the typical negative tail coming from the pion field (sea contribution). 
The proton form factor is in good agreement with the experiment whereas the 
neutron one is overestimated by a factor two. Both form factors show practically 
no dependence on the constituent mass M. For the proton the sea contribution 
is less than 5% which is not the case of the neutron one where the valence quarks 
and the sea contribute with similar magnitude and opposite signs. Concerning the 
neutron electric form factor one could expect that both the boson loop corrections 
as well as the incorporation of vector mesons 14 may modify it. As can be expected 
the isoscalar squared radius is reproduced well whereas the isovector one is by a 
factor of two too large (see Tab. 1). The Tab. 1 lists the Coulomb part of the 
proton-neutron mass splitting as well. 

In non-relativistic limit the axial form factor defined as usual by 



(N 2 (p) \Af(0)\ N^p')) = g A (q 2 ) - ^% <X 2 \a,r a \ X X ) (2.13) 



G E (q 2 ) = (N(p) \jr(0)\N(p')), 



(2.11) 



(2.12) 
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is also calculated 15 using the quantized nucleon state \Ni(p')). As can be seen 

from Fig. 5 the results are in good agreement with the experimantal dipole fit 

at finite momentum transfer Q 2 . The value at the origin (about 0.8), which is 

exactly the axial vector coupling constant, is lower than the experimental value 

1.2. However, it is a common problem of the effective chiral models. Kirchbach 

and Riska 16 relate it to the problem of a proper renormalization of the axial charge. 

Table 1 The nucleon squared radii computed for the constituent quark masses. Total, 
valence and sea contributions are given separately to see the influence of the 
vacuum polarization effects. Listed are also the N — A mass splitting as 
well as the Coulomb part of the neutron-proton mass splitting. In the last 
column the experimental values are given for comparison. 



Quantity 


Co 

320 MeV 


nstituent Quark M 

420 MeV 


ass 

465 MeV 


Experiment 


total 


sea 


total 


sea 


total 


sea 


< r 2 > T=0 Ifm 2 } 


0.66 


0.07 


0.61 


0.16 


0.53 


0.15 


0.62 


< r 2 > T=1 [fm 2 ] 


1.36 


0.71 


1.15 


0.74 


1.05 


0.73 


0.86 


< r 2 > n [fm 2 ] 


-0.35 


-0.32 


-0.27 


-0.29 


-0.26 


-0.29 


-0.12 


< r 2 > p [fm 2 ] 


1.01 


0.39 


0.88 


0.45 


0.79 


0.44 


0.74 


M A -M N [MeV] 


203 




261 




301 




294 


A EM [MeV] 


-0.87 


-0.18 


-1.00 


-0.73 


-1.09 


-0.33 


-0.76 ±0.30 



The sigma commutator calculated in the model comes out quite reasonable - 
about 35 MeV. 

We also studied possible modifications of the nucleon properties in medium 
(finite density) . The medium is simulated simply by varying the chemical potential 
/i in order to fix the medium density. The mass of the nucleonic soliton is depicted 
in Fig. 6 as function of the density. With increasing density the mass gets reduced 
as at some critical value the soliton disappears - a sign of a derealization of the 
nucleon. The mean radius shown in Fig. 6 supports this picture. At finite density 
the soliton gets extended (swelling) and less bound and close to some critical 
density it is getting delocalized. 
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3. NJL model with SU(3)-flavours 

The SU(3)-invariant NJ1 model withscalar and pseudoscalar coupling reads 

C N JL = q(x)W - m)q(x) - | [(q(x) X a q(x)) 2 + (q(x)i l5 X a q(x)) 2 ] (3.1) 

where m = diag(m u , m^, m s ) = mil + 7713 Ag is the quark mass matrix and A a 
are the usual Gell-Mann matrices with A = \J (§)1. Actually we use the current 
quark masses mg = \{tti u + m^). Integrating out the quarks the effective action 
is given by 

S eff = -Sp log (ip-m- g(a a X a + i-y 5 ir a \ a )) + ^-(a a a a + n a n a ) (3.2) 

Hereafter we work in the chiral limit and treat the strange quarkmass m s per- 
turbatively in first order. Similar to the SU(2) case we fix the parameters in 
the meson sector of the vacuum by reproducing the experimental values of the 
meson masses m n = !39MeV and mx = 496MeV, as well as the pion decay 
constant = 93MeV. This results in the relation m\lm\ = (m s + mo)/2mo 
with mo = \{jri u + m ( j). In order to fulfill it one needs a current mass mg of about 
6 MeV. To this end we use a proper time regularization with a modified weight 
function. The latter leads to a higher value of the E-commutator E = 51 MeV. 

For the SU(3) soliton we use the 5 , L r (2)-one Uq trivially embedded in SU(3) 

U,= { ;j (3.3) 

In order to quantize the classical soliton we follow the same scheme (cranking 
approxaimation) as in the SU(2) case (see eq.(2.8)) but with a SU(3) isorotating 
meson fields. As a next step we expand the effective action in fi^ up to second 
order 

L r0t = \l A B^B ~ ^|^8, (3-4) 

where Ij^q is the SU(3) tensor of the moment of inertia, as well as up to the first 
order in m s 

2 m s 



L^ = - 2 -^^e(i-D^(R)) (3.5) 
3m u + m d V 88 V V 
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and m s flA 

l(2) = _^ KabD W {r)SIb . (3 . 6) 

(8) 

All other corrections are nee; lected. Here D^(R) is the SU(3)-Wi erner function 
and by K^b we denote the tensor of the anomalous moments of inertia. These 
moments need no regularization, because they originate from the imaginary part 
of the effective Euclidean action and so they are finite. It should be stressed that 
the appearance of the linear term in fig in (3.4) is due to the discrete valence level 
in the spectrum whereas in the Skyrme model a term like this can be obtained 
only by adding the Wess-Zumino term. 

Now we quantize canonically by defining SU(3) right generators Ra=i,...,8 by 
Rj± = —dL/dVLj{, where L = L rot + iX 1 ) + L( 2 \ Using the collective hamiltonian, 
obtained in approximation 0(m%), in a strict perturbation theory to first order in 
m s we calculated the masses of the member of the octet and decouplet. We fixed 
m£* to the experimental value, because as in the SU(2)-case the total masses are 
always too large by a constant shift (zero-point energies are not subtracted). It 
should be stressed that we obtained the Gell-Mann Okubo relations 

2{m^ + ra^) = 3?tia + m S an d wiq — m^* = m^* — mj* = mj* — 
as well as the Guadignini formula 17 

m^* — m^* + m^v = (l/8)(H m A — 3 TO s) 

without any additional assumptions. The latter is derived in the Skyrme model 
only by adding a term associated to the hypercharge simply by hand. Our results 
are illustrated in Fig. 7 where the discrepancy between the theory and experiment 
for the masses of the members of the octet and decouplet is presented as a function 
of the strangness quark mass. 
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4. Summary 

The SU(2)-flavour Nambu - Jona-Lasinio model with scalar and pseudoscalar 
couplings and valence quarks combined with the semi-classical quantization (crank- 
ing approximation) provides quite a reasonable description of nucleon properties 
like masses, charge radii as well as electric and axial form factors. In the Sub- 
version the mass splittings between the octet and decuplet as well as within the 
multiplets are peproduced not only in a correct oder but also in a good agreement 
with the experimantal values. Thus, the assumed physical picture of the baryon as 
a soliton of three valence quarks coupled to the Dirac sea is supported. The results 
do not leave much room for additional meson degrees of freedom, in particular for 
vector mesons. Since the NJL model respects the chiral symmetry breaking but 
not confinement one can conclude that i) the chiral symmetry breaking mecha- 
nism is dominant for the nucleon structure and ii) the confinement seems to be 
not of great importance at least for the very low-energy properties considered. In 
particular, the (partial) restoration of the chiral symmetry occuring in medium 
modifies of the nucleon structure: the nucleon gets less bound and swelled and at 
some critical density it undergoes derealization. 

The work has been supported by the Bundesministerium fur Forschung und 
Technologie, Bonn, by the KFA Jillich (COSY - Project), the Deutsche Forschungge- 
meinschaft and by Polish Research Grant 2-0091-91-01. 
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FIGURE CAPTIONS 



1. The energy of the nucleon (thick full line), delta (thin full line) as well 
as the non-quantized hedgehog soliton energy (dashed line) as function of 
the constituent quark mass. 

2. The proton electric form factor as a function of the squared momentum 
transfer . The solid lines describe the form factor and the sea contribution 
(thin line) for the constituent quark mass 360 MeV. The dashed lines 
display the same information for the constituent quark mass 420 MeV. 

3. The neutron electric form factor as a function of the squared momentum 
transfer Q 2 . The solid lines describe the total form factor and the sea 
contribution (thin line) for the constituent quark mass 360 MeV. The 
dashed lines display the same information for the constituent quark mass 
420 MeV. 

4. The neutron charge distribution as a function of the radius. The solid 
bold line denotes the total charge, the dashed line denotes the valence 
contribution and the thin line stands for the sea contribution only. 

5. Soliton mass as a function of the baryon density for two different values 
of the constituent quark mass is compared with the mass of the soliton 
(dash-dotted line) of the medium modified cr-model. The separated con- 
tribuions coming from the valence quarks as well as from the polarization 
of the Fermi and Dirac sea for the constituent mass M = 420 MeV are 
also presented. 

6. Soliton mean square radius as a function of the baryon density for two 
different values of the constituent quark mass. 

7. The deviation of the theoretical masses from the experimental values as a 
function of the strangness quark mass for constituent quark mass M=420 
MeV. 
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